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Abstract 

Background: Duck enteritis virus (DEV) is an unassigned member in the family Herpesviridoe. To demonstrate 
further the evolutionary position of DEV in the family Herpesviridoe, we have described a 42,897-bp fragment. We 
demonstrated novel genomic organization at one end of the long (L) region and in the entire short (S) region in 
the Clone-03 strain of DEV. 

Results: A 42,897-bp fragment located downstream of the LOFR11 gene was amplified from the Clone-03 strain of 
DEV by using 'targeted gene walking PCR'. Twenty-two open reading frames (ORFs) were predicted and 
determined in the following order: S'-LORF11-RLORF1-ORF1-ICP4-S1-S2-US1-US10-SORF3-US2-MDV 
US5-US6-US7-US8-ORFx-US1-S2-S1-ICP4 -3'. This was different from that of the published VAC strain, both in the 
linkage of the L region and S region, and in the length of the US10 and US7 proteins. The MDV091.5-like gene, 
ORFx gene, SI gene and 52 gene were first observed in the DEV genome. The lengths of DEV US10 and US7 were 
determined to be 31 1 and 371 amino acids, respectively, in the Clone-03 strain of DEV, and these were different 
from those of other strains. The comparison of genomic organization in the fragment studied herein with those of 
other herpesviruses showed that DEV possesses some unique characteristics, such as the duplicated US1 at each 
end of the US region, and the US5, which showed no homology with those of other herpesviruses. In addition, the 
results of phylogenetic analysis of ORFs in the represented fragment indicated that DEV is closest to its 
counterparts VZV {Voricellovirus) and other avian herpesviruses. 

Conclusion: The molecular characteristics of the 42,897-bp fragment of Clone-03 have been found to be different 
from those of the VAC strain. The phylogenetic analysis of genes in this region showed that DEV should be a 
separate member of the subfamily Alphaherpesvirinoe. 
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Background 

Herpesviruses are among the most persistent of all 
pathogens because they have coevolved with their hosts 
over a long period of time, and they are relatively harm- 
less in immunocompetent hosts [1]. The family Herpes- 
viridae comprises approximately 100 members; these 
viruses infect a range of host species from humans and 
other mammals to birds, amphibians, and reptiles [2]. On 
the basis of differences in cellular tropism, genome orga- 
nization, and gene content, herpesviruses have been 
grouped into three subfamilies: Alphaherpesvirinae (a-), 
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Betaherpesvirinae (P-), and Gammaherpesvirinae (y-) 
[3,4]. Currently, duck enteritis virus (DEV), also known 
as duck plague virus (DPV) and duck herpesvirus- 1 [4], is 
an unassigned member of the family Herpesviridae [5]. 

Herpesviruses are enveloped viruses with a virion size 
over 100 nm [1]. The genomes of these viruses are lin- 
ear, double-stranded DNA, and they differ in size, 
sequence arrangements, and base composition [2]. They 
also vary significantly with respect to the presence and 
arrangement of inverted and directly repeated sequences 
[6]. Herpesvirus genomes differ in the arrangement of 
direct and inverted repeat regions with respect to 
unique regions. Six types of genome structures have 
been confirmed adequately in herpesviruses, which are 
designated by letters from A to F. The A type structure 



O© 201 1 Liu et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons 
BiolVlGCl CcntFc! I Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in 
any medium, provided the original work is properly cited. 



Liu et al. Virology Journal 201 1, 8:200 
http://www.virologyj.eom/content/8/1/200 



Page 2 of 1 2 



consists of a unique region flanked by a direct terminal 
repeat at the genome ends. Type B genomes contain 
variable numbers of a TR (terminal reiterations) at each 
end of the genome. In the C type genome, the number 
of direct terminal reiterations is small but sequences 
longer than 100 bp are directly repeated and subdivide 
the unique sequence of the genome into several well 
delineated stretches. The D type genome just has the 
repeated sequences at one terminus and in an inverted 
orientation internally. In the E group, the genome is 
divided into unique long (UL) and unique short (US) 
regions; each unique region is flanked by the inverted 
repeats. The sequences at the two termini of the F 
group are not identical and are not repeated directly or 
in an inverted orientation. It has been reported that 
DEV also contains linear, double-stranded DNA, and its 
genome was shown to be approximately 180 kb in size, 
with a G plus C content of 64.3% [7]. Genomic 
sequences of DEV have been reported recently by sev- 
eral Chinese research groups; however, discrepancies 
were found among these reports [8-18]. Genes in the 
UL region of DEV and their arrangement have been 
reported by our laboratory, and the results generally 
showed more similarity with Mardiviruses [8-13]. 
Another report showed that the LORF11 gene of the 
VAC strain is located at the leftmost end of the DEV 
genome, and that the LORF11 gene encoded a putative 
protein of 275 amino acids in the VAC strain [14]; both 
of these results differ from our previous results [12]. 
Meanwhile, several genes in the US region have also 
been reported [15-18]; however, the length of the puta- 
tive proteins encoded by the US 10 gene and US7 gene 
has been debated. In this study, we present a fragment 
of 42,897 bp, which contains one end of the L region 
that includes part of the LORF11 gene, which was 
absent from the published VAC strain, and the whole of 
the DEV S region. In addition, we demonstrated a dif- 
ferent genomic organization of the junction of the L 
region and the S region in this study. These results will 
provide a useful comparative dataset for the study of 
related genes in DEV and other herpesviruses. 

Results 

The features of the overall sequences and determination 
of ORFs 

A fragment of 42,897 bp downstream of the LORF11 
gene was amplified from the genome of the Clone-03 
strain of DEV in this study. The genome structure and 
the gene layout of this fragment are depicted in Figure 1. 
The fragment contained part of the sequence of the 
LORF11 gene [12], the rightmost part of the L region, the 
US region and its flanking sequences, and inverted 
repeats of the short region (IRS and TRS). The L region 
and IRS were interrupted by a set of tandem repeat 



sequences designated as a-type-like sequences [13], as in 
the case of the two regions in herpes simplex virus (HSV) 
[19]. Another a-type-like sequence was also found at the 
end of the TRS in the DEV genome. The overall G plus 
C ratio of the region sequenced was 46.09%. 

Twenty-two ORFs that contained more than 75 amino 
acids were found in the present study, which were in 
the order: 5'-LORFll-RLORFl-ORFl-lCP4-Sl-S2-USl- 
US1 0-SORF3- US2-MD V091 . 5- like- US3- US4- US 5- US6- 
US7-US8-ORFx-USl-S2-Sl-ICP4-3\ These ORFs were 
predicted to encode 17 putative proteins, with the 
exception of LORF11, because genes in the IRS and TRS 
were inverted and encoded the same proteins. The start 
locations of all ORFs were assumed to be the first possi- 
ble ATG. The motifs of each ORF are listed in Table 1. 

The confirmation of the junction between the L region 
and the S region 

Owing to the different linkages of the L region and S 
region found in the genome sequences of the published 
DEV VAC strain [14] and our above-described sequence 
in the Clone-03 strain of DEV, a pair of specific primers 
was designed to confirm the junction of the L region 
and S region in the DEV genome. The forward primer, 
L25, was located in the LORF11 gene (GenBank no. 
EU294364), which is a gene in the DEV UL region that 
had only one copy in the genome compared with the 
genomes of other alphaherpesviruses. The reverse pri- 
mer, L26, was located in the SORF3 gene, which is a 
gene in the US region of the DEV genome that also has 
a single copy in the DEV genome. The PCR product 
was used as the model for the second nested PCR after 
dilution to 1 in 1,000. We obtained four different frag- 
ments (Figure 1), and they were 4,553 bp, 4,689 bp, 
4,743 bp, and 5,547 bp in length, respectively. The 
results of sequencing of the four fragments showed that 
they were parallel with the sequences obtained using 
'targeted gene walking PCR'. Consequently, we deter- 
mined that the linkage between the L region and the S 
region should be in the following order: S'-LORFll- 
RLORF1-ORF1-ICP4 -Sl-S2-USl-US10-SORF3-3\ 

A 207-bp insertion in both the IRS and the TRS regions 
was not found in their counterparts in the DEV VAC 
strain 

In addition to the linkage of the L region and the S 
region, two insertions of 207 bp were found in the pre- 
sented fragment in both the IRS region and the TRS 
region (Figure 2), when compared with the published 
VAC genome. The 158-bp sequence at the 3' end of the 
207-bp sequence of the IRS region was complemented 
with a fragment of the same length at the 5' end of the 
207-bp sequence of the TRS region. The remaining 
49-bp fragment in each of the insertions was dissociated 
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Clone-03 strain 




42,897bp 



Figure 1 The comparison of genomic organization between the DEV Clone-03 strain and the published VAC strain, together with the 
PCR strategy. The upper part shows the genomic organization of sequences corresponding to those in the present study in the published 
sequences of the genome of the DEV VAC strain. The genomic organization of the presented fragment in the DEV Clone-03 is listed below. The 
red and blank arrows indicate the ORFs, with different colours to make two adjacent ORFs evident. The pink arrows indicate the ORFs that were 
not detected in other heprsviruses but were detected in DEV. The dark blue indicates the UL region; the dashed gray boxes in the genome of 
DEV Clone-03 indicate uncertain regions. The PCR strategy used to obtain the fragment is depicted at the bottom. The blue arrows indicate the 
twelve overlapping fragments that were obtained to form a continuous DNA fragment; the primer positions and amplification directions are also 
shown. The arrows indicate the genes and the interval shows the relative position of the two adjacent ORFs. The confirmation PCR strategy is 
depicted in the centre. The bars indicate the PCR product with the primers embedded. The green ellipse indicates the predicted origins of 
replication (oriS). 

V. J 



and not complemented. Both of the fragments were rich 
in A plus T, with a content of 67.15%. 

The characteristics of new ORFs detected in the fragment 

Two ORFs, designated RLORF1 and ORF1, were 
detected in the region upstream of the S region. 
Another copy of ORF1 was found to the left of the DEV 
L region of the genome [13]. The RLORF1 and ORF1 
encoded two putative proteins of 109 and 81 amino 
acids, respectively. Four phosphorylation sites were pre- 
dicted in the sequence of RLORF1. 



In addition, eight ORFs encoding four different puta- 
tive proteins (SI, S2, ICP4, and US1) in the RS region 
were detected. Of these proteins, SI and S2 were identi- 
fied for the first time in the present study. The SI gene 
encoded a putative protein of 92 amino acids, and four 
phosphorylation sites were predicted. No homologue of 
SI was found in the proteins encoded by other herpes- 
viruses. Another unique gene in the RS region was S2, 
which encoded a putative protein of 96 amino acids that 
contained just six phosphorylation sites. ICP4 and US1 
were the same as previously described [13,14]. 
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Table 1 Core promoters searched in the neural network and polyadenylation signals predicted by POLYADQ 



Gene 


Promoter location 3 


Promoter score 


TATA sequence 


TATA location 


Poly(A) sequence 


Poly(A) location 


Poly(A) score 


RL0RF1 


NP b 


NP 


NP 


NP 


NP 


NP 


NP 


0RF1 


5,564-5,613 


0.93 


ATATAAAGCGGTAGT 


5,575-5,589 


NP 


NP 


NP 


ICP4 


1 0,928-1 0,977r c 


0.88 


^GTAAAAT 


1 0,960-1 0,969r 


AATAAA 


5,867-5,872r 


0.317475 


SI 


13,570-1 3,61 9r 


0.85 


CTATCTAAGGCGACC 


13,602-13,61 1r 


NP 


NP 


NP 


S2 


14,725-14,774 


1.00 


NP 


NP 


NP 


NP 


NP 


US1 


15,606-15,655 


0.99 


GCCTAAAAAGCACCG 


15,613-15,628 


AATAAA 


17,015-17,020 


0.644949 


US10 


1 7,003-1 7,052r 


0.94 


CAATAAACACCGOT 


17,014-17,028 


NP 


NP 


NP 


S0RF3 


1 9,280-1 9,329r 


0.99 


GC^AAAAG 


19,31 3-1 9,322r 


AATAAAr 


1 8,248-1 8,253r 


0.138679 


US2 


20,407-20,456r 


0.85 


GTCTAAAAGGCAGAG 


20,434-20,448r 


NP 


NP 


NP 


MDV091.5-like 


NP 


NP 


NP 


NP 


NP 


NP 


NP 


US3 


20,299-20,348 


0.85 


CCCATAAATG 


20,305-20,314 


NP 


NP 


NP 


US4 


21,799-21,848 


0.96 


GTATAAATOGACAA 


21,807-21,821 


AATAAA 


23,310-23,315 


0.385382 


US5 


Zj,jJJ Zj,jO^ 


0 88 


vJlv-l 1 VJ 1 VJ 1 1 1 rA 1 rA 1 


Zj,jjj ZjjjM-? 


AATAAA 

/AAA I AAAA/A 


Z J, I / U Z J, I / J 


u.zuu^o/ 


US6 


25,058-25,107 


0.94 


CGGCAATATGTATAT 


25,062-25,076 


NP 


NP 


NP 


US7 


26,158-26,207 


0.92 


ATATAATOCTACGC 


26,167-26,181 


NP 


NP 


NP 


US8 


NP 


NP 


NP 


NP 


NP 


NP 


NP 


ORFx 


29,414-29,463 


0.93 


GTATATOGGCCGAC 


29,421-29,435 


NP 


NP 


NP 


US1 


3 1,386-31, 435 r 


0.99 


GCCTAAAAAGCACCG 


31,41 3-3 1,427r 


AATAAA 


30,021-30,026 


0.644949 


S2 


32,267-32,3 16r 


1.00 


NP 


NP 


NP 


NP 


NP 


SI 


33,422-33,471 


0.85 


CTATCTAAGGCGACC 


33,431-33,445 


NP 


NP 


NP 


ICP4 


36,064-36,113 


0.88 


mGTAAAAT 


36,072-36,081 


AATAAA 


41,170-41,175 


0.317475 



a the position is according to the sequence of the whole fragment of 42,897 bp. 
b NP indicates no prediction. 
c r indicates reverse direction. 



r \ 
RLORF1 MDV091.5-like 




5 s - 

T AT AT GC T AAC AT T T AC C C G C C C T AAGAAAA 
CTTACCCATAACCACCGATAATGGTATTCAT 
TCAGGAAATAGTAAGAATAACATTATTGGAT 
ATGGTTTAATTATATGCTAACATTTACCCGC 
C C T AAGAAAAC T T AC C C AT AAAC AC T GT T AA 
TTATATTCAACCGGGAACTAGTAGAAATAAT 
ATTACTGGATCTAGCCGTAAT-3 N 

Figure 2 The position of the two 207-bp insertions in the IRS and TRS regions in DEV Clone-03 strain comparing to the VAC strain. 

The red and blank arrows indicate the main genes in the IRS and TRS regions. The gray boxes indicate the interval between the ICP4 and US1 
genes. The blue and green boxes with dashed lines indicate the 207-bp fragment that was deleted from the IRS and TRS regions in the 
published VAC strain, respectively. Letters a and b indicate the fragment in the IRS and TRS regions, respectively. The sequences of the two 
fragments are listed in the boxes at the bottom. 



5^- 

TAATTAACAGTGTTTATGGGTAAGTTTTCTT 
AGGGCGGGTAAATGTTAGCATATAATTAAAC 
CATATCCAATAATGTTATTCTTACTATTTCC 
TGAATGAATACCATTATCGGTGGTTATGGGT 
AAGTTTTCTTAGGGCGGGTAAATGTTAGCAT 
ATAATTACGGCTAGATCCAGTAATATTATTT 
CTACTAGTTCCCGGTTGAATA-3 N 
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The DEV US region contained 11 ORFs that were likely 
to code for 11 proteins (Figure 1), which included homolo- 
gues of the HSV-1 genes US10, US2, US3, US4, US6, US7 
and US8 [20]. Interestingly, a unique ORF in the DEV US 
region, located downstream of US8, was predicted in the 
present study and named ORFx. The ORFx encoded a 
putative peptide of 118 amino acids. One transmembrane 
domain was detected in the ORFx between residue posi- 
tions 95 and 115 at the N-terminus. Remarkably, the 
length of our DEV US 10 was 311 amino acids, which was 
different from published results of 168, 169 and 298 
amino acids [14,15,17]. We also found a sequence of 13 
amino acids, CSFWCCLGHAATC (Additional file 1, Fig- 
ure SI), which mapped to amino acids 236-248 and con- 
formed to the C-C-H-C zinc finger motif as described in 
equine herpesvirus-1 (EHV-1) [21,22] 

A new gene was predicted in this study, which was 
327 bp in length and overlapped 197 bp at the 3'-termi- 
nus of the US2 gene. It was homologous to the proteins 
encoded by Marek's disease virus-1 (MDV-1), MDV-2 
and HVT and was designated MDV091.5-like gene. 
BLAST searches using the amino acid sequence showed 
that this protein had some amino acid similarity with 
putative nucleotide-binding oligomerization domain- 
containing protein 2 of Gasterostrus aculeatus, the puta- 
tive lyase of Rhodococcus erythropolis, and bacterial 
valyl-tRNA synthetase. 

The transmembrane regions of the proteins encoded 
by the genes in the presented fragment are depicted in 



Figure 3. The conserved domains of US1, SORF3, US2, 
US3, US4, US6, US7, US8 proteins are shown in Addi- 
tional file 2, Figure S2, Additional file 3, Figure S3, 
Additional file 4, Figure S4, Additional file 5, Figure S5, 
Additional file 6, Figure S6, Additional file 7, Figure S7, 
Additional file 8, Figure S8, Additional file 9, Figure S9, 
Additional file 10, Figure S10, respectively. 

Phylogenetic analysis 

Phylogenetic rooted trees were constructed from align- 
ments of the putative proteins with their homologues in 
other alphaherpesviruses and are shown in Figure 4 and 
5. The DEV US2 gene, US3 gene, US6 gene, US7 gene 
and US 10 gene showed closer relationships with mem- 
bers of Mardivirus. However, US1 showed a closer rela- 
tionship between DEV and members of Simplexvirus 
and Varicellovirus. The DEV US4 gene showed more 
similarity with infectious laryngotracheitis virus (ILTV), 
and both clustered into the subfamily Varicellovirus 
(Figure 4). The DEV US8 gene fell into an outgroup 
position with respect to members of subfamily Alpha- 
herpesvirinae (Figure 5), which implies that a recombi- 
nation event may have occurred during the origin and 
evolution of the virus. 

The comparison of gene layouts in the US region of DEV 
with those in other alphaherpesviruses 

A comparison of the genetic organization of selected 
alphaherpesvirus US segment genes is presented in 



SORF3 

Cytoplasmic side C-terminus 



Outer membrane side <: — N-terminus 




US4 



Cytoplasmic side 

C-terminus 

N-terminus 



Outer membrane side 



US5 



Cytoplasmic side 




N-terminus 
Outer membrane side 



US6 



Cytopasmic side 
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Outer membrane side 
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Outer membrane side 
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Figure 3 Transmembrane regions in the proteins encoded by DEV genes in the present fragment. The TMpred server was used to 
predict transmembrane regions. The dashed purple lines indicate the amino acid sequence of a potential signal peptide that may be cleaved 
during the maturation of proteins. 
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Figure 4 Evolutionary relationship of US1 7 US2 7 US3 and US4 proteins in DEV S regions within the subfamily Alphaherpesvirinae. The 

phylogenetic tree was generated using the MEGALIGN (DNAStar) program. Sequence distances indicated by the scale were calculated using the 
Gonnet 250 matrix in LASERGENE. 



Figure 6. Despite obvious similarities, there were marked 
differences in gene content, organization and localiza- 
tion between DEV and other alphaherpesviruses. Never- 
theless, these overall gene layouts are consistent with a 
model that accounts for the divergence of alphaherpes- 
virus from a common ancestor by a number of homolo- 
gous and semihomologous recombination events, which 
resulted in concomitant loss or gain of US genes [23]. 

Origins of replication in the S region 

Two well-defined origins of replication were found in 
the IRS and TRS of the DEV genome, designated oriS. 
The two oriS were palindromic structures and contained 
the same sequence features: two inverted 9-bp 
sequences, which were identical to that recognized by 
the origin-binding protein (OBP) encoded by the UL9 
binding sequence (GTTCGCAC), separated by a 43-bp 
AT-rich spacer sequence (76.75% A+T) (Figure 7). The 
features were the same as described for PRV (Pseudora- 
bies virus) [24] and equine herpesvirus- 1 (EHV-1) [25]. 



Discussion 

Our laboratory has been engaged for many years in ana- 
lyzing the genome sequences of DEV [8-13]. After we 
had completed the genome sequence of DEV Clone-03, 
a DEV VAC genome sequence was also published by 
other researchers [14]. However, some differences were 
detected by comparison of parts of our DEV Clone-03 
strain with those of the DEV VAC strain. Herein, we 
presented the sequence of a 42,897-bp fragment 
anchored in the LORF11 gene of the DEV genome 
which was located at the rightward end of the UL region 
[12], by using the method of 'targeted gene walking 
PCR' (Figure 1). Comparison of the sequence of the 
fragment with that of the DEV VAC strain showed that 
our Clone-03 strain of DEV had a different gene order 
from that of the DEV VAC strain in this region. Conse- 
quently, we designed an additional four pairs of primers 
according to the new sequences and confirmed the 
result using nested PCR (Figure 1). The two methods 
obtained the same sequences, and it was demonstrated 
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Figure 5 Evolutionary relationship of US6, US7, US8 and US10 proteins in DEV S regions within the subfamily Alphaherpesvirinae. The 

phylogenetic tree was generated using the MEGALIGN (DNAStar) program. Sequence distances indicated by the scale were calculated using the 
Gonnet 250 matrix in LASERGENE. 



that the genes in this region should be in the following 
order: 5'-LORFll-RLORFl-ORFl-ICP4-Sl-S2-USl-US10- 
SORF3-3', which is different from the DEV VAC strain, 
in which the gene order is 5'-LORFll-UL-ICP4-USl- 
US10-SORF3-3' [14]. The different linkage pattern 
between DEV Clone-03 and the VAC strain in the L 
region and S region is difficult to explain and requires 
further investigation, although a different linkage 
between the L and S regions of HSV was also observed 
between wild-type virus and cell-adapted virus [26,27]. 

Interestingly, we also found some novel characteristics 
of the sequences in the S region of the Clone-03 strain 
of DEV. Two insertions of 207 bp in the IRS and TRS 
regions were found in the DEV Clone-03 strain that 
were absent from the VAC strain. It has been reported 
that some fragments were lost during serial passage of 
MDV [28]. Hence, we speculated that the insertion of 
the two 207-bp fragments in the DEV Clone-03 strain 
and their absence from the VAC strain might be due to 



the different passage levels [28]. The SI gene, S2 gene, 
RLORF1 gene, ORF1 gene and ORFx gene that were 
observed in the Clone-03 strain in this study also had 
similar sequences in the VAC genome; however, those 
genes showed no homologues in other alphaherpes- 
viruses. Those genes may be potential markers to differ- 
entiate DEV from other alphaherpesviruses. 

Davison and McGeoch concluded that differences in 
gene layout in the S component between HSV-1 and 
VZV have resulted from expansion and contraction of 
IRS/TRS during evolution [23]. This may also be the case 
for the DEV genome. Unlike those of MDV-1, MDV-2 
and HVT, the DEV US1 gene was duplicated and also 
inverted to the other end of the US, as is that of PRV 
[24]. Similarly, the presence of two copies of the US1 
gene in DEV does not imply that the virus expresses two 
forms of ICP22 [24]. Although the pattern of the two 
copies of the US1 gene in the DEV genome showed a 
similar gene layout to those of PRV, the existence of the 
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Figure 6 Comparison of genes in the DEV and alphaherpesvirus US regions. The comparison was based on the published sequences (37, 
21, 20, 49). The red and blank arrows indicate the ORFs that were homologous in different alphaherpesviruses. The gray arrows indicate the 
ORFs that not exist in all the selected viruses. The green bars in the leftmost and rightmost parts indicate the sequences of the genomes, except 
for the US region. The targeted regions compared in this study are boxed. 



LORF11 gene at the rightward end of the UL region indi- 
cated that the organization of the DEV genome may be 
similar to that of other avian herpesviruses. The presence 
of the SORF3 gene and the MDV09L5-like gene, and the 
translocation of the US10 gene in the DEV genome, 
further suggests a close relationship between DEV and 
other avian herpesviruses. In addition, the phylogenetic 
analysis of most genes in the presented fragment further 
indicated a close relationship between DEV and viruses 
in the subfamily Mardivirus. However, the US region of 
DEV contained some genes that were absent from the 
genomes of other avian herpesviruses, such as US4 and 
ORFx, which indicates that DEV may be a unique mem- 
ber of the subfamily Alphaherpesvirinae. 

Replication of the viral genome is a central event in the 
life cycle of herpesviruses. The initiation of viral DNA 
synthesis marks the commitment of the infected cell to 
the production of new infectious virus and, in most 
instances, cell death. HSV-1 contains three origins of 
DNA replication of two types: one copy of oriL located at 
the centre of the UL region of the genome and two copies 
of oriS located in the repeat regions that flank the US 
region of the genome [29]. The reasons for the three 
potential origins of replication in the viral genome are not 
apparent in HSV. In this study, we predicted two copies of 



oriS in the RS region of DEV. It has been reported that the 
deletion of the oriL in HSV resulted in reduced replication 
in mouse tissues and reduced reactivation from latent 
infection. Thus, oriL may be required for DNA replication 
in certain tissues [29] . Although oriL was absent from the 
DEV genome, the core sequence of oriS, which typically 
contains an origin recognition element and a DNA- 
unwinding element, was unchanged [29]. This absence of 
oriL from DEV may be associated with the evolution of 
the viral genome, may lead to different characteristics of 
the replication of DEV from those of other herpesviruses, 
and may even result in functional deletions from the gen- 
ome of DEV in comparison with other herpesviruses. 

Conclusion 

In this study, we demonstrated a different organization 
of genes in the rightward part of the L region and the 
whole S region in the Clone-03 strain of DEV, when 
compared with the VAC strain. Several novel character- 
istics were also detected in this region that have not 
been reported in the VAC strain, including the presence 
of SI, S2, ORFx and MDV091.5-like genes and two 
insertions in the IRS and TRS regions. The genomic 
order and the characteristics of the genes in this region, 
together with phylogenetic analysis based on the 
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putative proteins encoded by the genes investigated in 
the present study showed that DEV should be a unique 
member of the subfamily Alphaherpesvirinae. 

Methods 

Virus stock preparation 

The Clone-03 strain of DEV was used in this study 
[8-13]. The virus stocks were produced by propagation 
in chicken embryo fibroblasts (CEF) in Dulbecco's mini- 
mum essential medium (DMEM) with 8% fetal bovine 
serum. The infected CEFs were harvested when the 
cytopathic effect (CPE) reached 80%. After three freeze- 
thaw cycles, the virus stocks were confirmed primarily 
by electron microscopy and polymerase chain reaction 
(PCR) as described previously [8-13]. 

DNA extraction, polymerase chain reaction and 
sequencing 

The viral DNA was extracted from the virus stocks as 
described previously [9] . The 'targeted gene walking PCR', 
as described previously [30,3 1], was used to amplify the 
targeted DEV genome fragment, as illustrated in Figure 1. 
Briefly, four nonspecific walking' primers, Nl, N2, N3 and 
N4 [10] were used to walk the genome of DEV. A pair of 
specific primers, LI and L2, was designed on the basis of 
the partial sequence of DEV LORF11 published in Gen- 
Bank (GenBank no. EU294364) [12]. The PCR was carried 
out by using LI and L2 as forward primers; the four non- 
specific primers were used as reverse primers. Finally, a 
2,998-bp fragment (Fl) was amplified, anchored from the 
DEV LORF11 gene. Targeted primer L3 and internal 



primer L4 were designed on the basis of the newly gener- 
ated fragment, and was used to amplify the neighbouring 
gene fragment with one of the four nonspecific primers. 
Similarly, primers L5-L6, L7-L8, L9-L10, L11-L12, L13- 
L14, L15-L16, L17-L18, L19-L20, L21-L22 and L23-L24 
(Figure 1) were designed and used in the subsequent PCR 
amplifications. The primers used in PCR amplifications in 
this study are listed in Table 2. 

The PCR was carried out in a 25 ul reaction volume 
as described previously [9]. The reaction was performed 
at 95°C for 5 min, followed by 30 cycles of 94°C for 1 
min, 50°C for 1 min and 72°C for 3 min; the reaction 
was ended by elongation at 72°C for 10 min. The PCR 
products were analyzed on a 0.8% agarose gel. The PCR 
products were sequenced directly or cloned into the 
pMD18-T vector (TaKaRa, Dalian, China) according to 
the manufacturers instructions and used for sequencing. 
Each of the fragments was sequenced at least three 
times from different PCR products. 

The determination of open reading frames (ORFs) in the 
presented fragment and genomic organization in the 
junction between the L region and S region in the DEV 
genome 

The sequences obtained were assembled using the Gene 
Runner (version 3.00, Hastings Software, Inc., Hudson, 
NY, USA). The ORFs and genomic organization in the 
junction of the L region and S region, and the layout of 
genes in the S region, were determined by comparison 
with the sequence counterparts of Marek's disease virus 
(MDV), HSV-1 and varicella-zoster virus (VZV). The 
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Table 2 Sequences of oligonucleotides primers used for PCR amplification 


Primer 


Sequences (5'-3') 


Position 3 


L1 


AGTCCAGTCATCTCCATCCG 


1,562-1,581 


L2 


ACGA^GGCTGTGCTGTAG 


1,730-1,749 


L3 


CTGTC^AAGG^AGGGCTGGC 


4,532-4,553 


L4 


GTAGGAAATA^GAGCCGAG 


4 826-4 845 


L5 


TTCTGATG I I I I GGCAAGCC 


7,578-7,597 


L6 


CAGAAJGGCGCYTTGTTTGG 


7,669-7,688 


L7 


^GAAGATAGG^GCTCGTAG 


1 1,272-1 1,292 


L8 


ATACAGGAAAA^AACGAT 


1 1 367-1 1 385 


L9 


ATGTAGCAG^G^CAAAC 


1 4,1 93-1 4,21 2 


L10 


AAGATGAGTCAACACCGAAGG 


14 ; 445-14 ; 465 


L1 1 


TA^CCATCCAG^GCTCCC 


1 7,732-1 7,751 


L12 


C^GTAAAGCTGGCCGCTAC 


1 7,978-1 7,997 


L13 


CGATCTGC^CGC^CCG 


21,725-21,744 


L14 


TAGCTGGTATGGCAACAATG 


21 895-21 914 


L15 


GAAG^AACGGAGGAAGTA^ 


26 043-26063 


L16 


GCATATAA^ACTACGCAACC 


26,165-26,185 


L17 


GTCATCCTTGTTATGTTGA 


29,697-29,715 


L18 


CCTAC^GGTGGTCGGCCA 


2Q 795-29 81 3 


L19 


CAGGA^GATAACTAACC 


34 266-34 284 


L20 


ATAAGCGCACTAGATGGCAG 


34 473-34 492 


L21 


TGAACGGACC^GCTAATGAC 


37,793-37,813 


L22 


GAGGGGTGGTACTGG^CCG 


38,120-38,139 


L23 


CCTACAATAACCTGGGAACT 


40,770-40,789 


L24 


GATC^GTCCGATGGGGATG 


40 949-40 968 


L25 


ATGGGACAGTCCCTACCG^GGCCTCGA^CAAAGCTOTCAG 


1-43 


L26 


GGACTGCAGGCC^TCAACGGCCCCTGATATGGCTACTATGTC 

vJ vJ/\v_ i vj v.r\vjvjv„v. i i i v^nnv_vjvjv.v_v.v. i vj/ \ i / \ i vj vj v_ i /\v_ i n i vj i v_. 


1 8 716-1 8 758r b 


L27 


GCCAGCCCTAACCTOAGACAG 


4,532-4,553r 


L28 


CAACCCCGCCCACAATAAGC 


4,043-4,062 


L29 


GCTOCATGC^CCCCGC 


8,71 2-8,731 r 


L30 


CGAACCGTCACAGTCTGCAG 


8,592-8,61 1 


L31 


GCCGCGCAGTAAGTACTCAG 


13,315-13,334 


L32 


TGCCGATATCATOGTOAT 


13,212-13,231 


N1 


TATAGGm(C/A)TG^ 


NP C 


N2 


C^GGAGCTG 


NP 


N3 


GAATGTGA(A/G)AA 


NP 


N4 


CATGTCTGCCGA 


NP 



a the position is according to the sequence of the whole fragment of 42,897 bp. 
b r indicates the reverse direction. 
c NP indicates no prediction. 



same program was used to detect ORFs encoding pro- 
teins of greater than or equal to 75 amino acids with a 
methionine (M) start codon. The predicted ORFs and 
flanking sequences were evaluated for coding potential 
by detecting the promoter http://www.fruitfly.org/seq_- 
tools/promoter.html[32], and the presence of TATA box 
http://motif.genome.jp/ and transcription terminal sig- 
nals http://rulai.cshl.org/tools/polyadq/polyadq_form. 
html. Searches of the deduced proteins for signal pep- 
tides http://www.cbs.dtu.dk/services/SignalP/, transmem- 
brane regions http://www.ch.embnet.org/software/ 



TMPRED_form.html, N-linked glycosylation sites http:// 
www.cbs.dtu.dk/services/NetNGlyc/ and serine, threo- 
nine and tyrosine phosphorylation sites http://www.cbs. 
dtu.dk/services/NetPhos/ were also performed online. 
The secondary structure of sequences in the oriS was 
constructed by using GeneQuest in DNAStar. 

Confirmation of the junction between the L region and 
the S region by specific PCR 

Owing to the different order of genes in the junction of 
the L region and S region in the DEV Clone-03 in this 
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study and the reported DEV VAC strain [14], one pair 
of specific primers, L25 and L26 (Table 2), was designed 
to confirm the result. Primer L25 was located within the 
LORF11 gene and L26 was located within the SORF3 
gene. This pair of primers was used in the first nested 
PCR. Other primers, L27, L28-L29, L30-L31, and L32 
(Table 2), were also used in the second nested PCR. 
The position of the primers and the strategy for confir- 
mation of the sequence are shown in Figure 1. 

The PCR was carried out in a 25 (il reaction volume. 
The first nest of the PCR reaction was performed at 95° 
C for 5 min, followed by 35 cycles of 94°C for 1 min, 
50°C for 1 min and 72°C for 8 min; the reaction was 
ended by elongation at 72°C for 10 min. The PCR pro- 
duct was analyzed on a 0.8% agarose gel and was used 
as the template for the second nest. The second nested 
PCR was performed at 95°C for 5 min, followed by 30 
cycles of 94°C for 1 min, 53°C for 1 min and 72°C for 3 
min; the reaction was ended by elongation at 72°C for 
10 min. The products of the second nested PCR were 
cloned and sequenced, respectively. 

Phylogenetic analysis 

Homologue searches were conducted using BLAST 
searching [33] and phylogenetic analysis was performed 
using the MEGALIGN program in Lasergene (DNAStar) 
with CLUSTAL W multiple alignment and weight 
matrix Gonnet 250 [13]. The result was confirmed by 
use of the MAGE package (Version 4.0). The sequences 
of the herpesviruses that were used as reference strains 
for homology analysis were obtained from the GenBank 
database and the GenBank accession numbers are given 
in the phylogenetic trees. 

GenBank accession numbers 

The DNA sequence of 42,897 bp from the DEV Clone- 
03 genome has been deposited in the GenBank database 
with the GenBank accession no. HQ009801. 

Additional material 



Additional file 5: Figure S5: Multiple alignments of homologues 
based on the amino acid sequences in the N-terminus of US3 
proteins of DEV Clone-03 and other typical strains of the subfamily 

Alphaherpesvirinae. The conserved domains (from I to VI) were 
indicated by pink boxes. 

Additional file 6: Figure S6: Multiple alignments of homologues 
based on the amino acid sequences in the C-terminus of US3 
proteins of DEV Clone-03 and other typical strains of the subfamily 

Alphaherpesvirinae. The conserved domains (from VII to XI) were 
indicated by pink boxes. 

Additional file 7: Figure S7: Multiple alignments of homologues 
based on US4 proteins of DEV Clone-03 and other typical strains of 
the subfamily Alphaherpesvirinae. The conserved domains were 
indicated by pink boxes. 

Additional file 8: Figure S8: Multiple alignments of homologues 
based on US6 proteins of DEV Clone-03 and other typical strains of 
the subfamily Alphaherpesvirinae. The conserved domains were 
indicated by pink boxes. 

Additional file 9: Figure S9: Multiple alignments of homologues 
based on US7 proteins of DEV Clone-03 and other typical strains of 
the subfamily Alphaherpesvirinae. 

Additional file 10: Figure S10: Multiple alignments of homologues 
based on US8 proteins of DEV Clone-03 and other typical strains of 
the subfamily Alphaherpesvirinae. The conserved domains were 
indicated by pink boxes. 



List of abbreviations used 

DEV: duck enteritis virus; DVE: duck viral enteritis; DPV: duck plaque virus; 
ORF: open reading frame; L: long; S: short; UL: unique long; US: unique short; 
a: Alphaherpesvirinae; (3: Betaherpesvirinae; y. Gammaherpesvirinae; CEF: 
chicken embryo fibroblasts; DMEM: Dulbecco's minimum essential medium; 
CPE: cytopathic effect; PCR: polymerase chain reaction; MDV: marek's disease 
virus; HVT: turkey herpesvirus; VZV: varicella-zoster virus; PRV: pseudorabies 
virus; M: methionine; ILTV: infectious laryngotracheitis; EHV: equine 
herpesvirus; CeHV: cercopithecine herpesvirus; CaHV-1: canid herpesvirus-1. 

Authors' contributions 

XL, SL and XK designed research; XL, ZH, YS and YL performed research; XL, 
SL, XK and HL analyzed data; and XL, SL and XK wrote the paper. All authors 
read and approved the final manuscript. 

Competing interests 

The authors declare that they have no competing interests. 

Received: 14 October 2010 Accepted: 2 May 201 1 
Published: 2 May 201 1 

References 

1. E Gershburg, JS Pagano, Conserved herpesvirus protein kinases. Biochim 
Biophys Acta. 1784, 203-212 (2008) 

2. JP Weir, Genomic organization and evolution of the human herpesviruses. 
Virus Genes. 16, 85-93 (1998). doi:1 0.1 023/A:1 00790591 0939 

3. B Roizman, LE Carmichael, F Deinhardt, G De-The, AJ Nahmias, W Plowright, 
F Rapp, P Sheldrick, M Takahashi, K Wolf, Herpesviridae. Definition, 
provisional nomenclature, and taxonomy. The herpesvirus study group, the 
international committee on taxonomy of viruses. Intervirology. 16, 201-217 
(1981). doi:1 0.1 159/000149269 

4. B Roizmann, RC Desrosiers, B Fleckenstein, C Lopez, AC Minson, MJ 
Studdert, The family Herpesviridae: an update. Arch Virol. 123, 425-449 
(1992). doi:10.1007/BF01317276 

5. AJ Davison, R Eberle, B Ehlers, GS Hayward, DJ McGeoch, AC Minson, PE 
Pellett, B Roizman, MJ Studdert, E Thiry, The order Herpesvirales. Arch Virol. 
154, 171-177 (2009). doi:1 0.1 007/s00705-008-0278-4 

6. B Roizman, The human herpesviruses in The herpesviridae, a brief 
introduction, in Fields of Virology, ed. by Roizman B, Whitley RJ, Lopez C 
(New York: Raven Press, 1993), pp. 1-9 



Additional file 1: Figure SI: Multiple alignments of homologues 
based on US10 proteins of DEV Clone-03 and other typical strains 
of the subfamily Alphaherpesvirinae. The pink box indicate the 
probable C-C-H-C zinc finger motif in US10 proteins by comparison with 
their homologues in other herpesviruses. 

Additional file 2: Figure S2: Multiple alignments of homologues 
based on US1 proteins of DEV Clone-03 and other typical strains of 
the subfamily Alphaherpesvirinae. 

Additional file 3: Figure S3: Multiple alignments of homologues 
based on SORF3 proteins of DEV Clone-03 and other avian 
herpesviruses. 

Additional file 4: Figure S4: Multiple alignments of homologues 
based on US2 proteins of DEV Clone-03 and other typical strains of 
the subfamily Alphaherpesvirinae. The conserved domains were 
indicated by pink boxes. 



Liu et al. Virology Journal 201 1, 8:200 Page 1 2 of 1 2 

http://www.virologyj.eom/content/8/1/200 



10. 



12. 



17. 



19. 



20. 



R Gardner, J Wilkerson, JC Johnson, Molecular characterization of the DNA 
of anatid herpesvirus 1. Intervirology. 36, 99-112 (1993) 
H Li, S Liu, Z Han, Y Shao, S Chen, X Kong, Comparative analysis of the 
genes UL1 through UL7 of the duck enteritis virus and other herpesviruses 
of the subfamily alphaherpesvirinae. Genet Mo I Bio. 32, 121-128 (2009). 
doi:1 0.1 590/S1 41 5-47572009005000003 

H Li, S Liu, X Kong, Characterization of the genes encoding UL24, TK and 
gH proteins from duck enteritis virus (DEV): a proof for the classification of 
DEV. Virus Genes. 33, 221-227 (2006). doi:1 0.1 007/s1 1262-005-0060-6 
S Liu, S Chen, H Li, X Kong, Molecular characterization of the herpes 
simplex virus 1 (HSV-1) homologues, UL25 to UL30, in duck enteritis virus 
(DEV). Gene. 401, 88-96 (2007). doi:10.1016/j.gene.2007.06.022 
R An, H Li, Z Han, Y Shao, S Liu, X Kong, The UL31 to UL35 gene 
sequences of duck enteritis virus correspond to their homologs in herpes 
simplex virus 1. Acta Virol. 52, 23-30 (2008) 

SW Liu, HX Li, Y Li, ZX Han, YH Shao, R An, XG Kong, Pylogeny of duck 
enteritis virus: Evolutionary relationship in the Family Herpesviridae. 
Intervirology. 51, 151-165 (2008). doi:1 0.1 159/000146865 
XL Liu, SW Liu, HX Li, ZX Han, YH Shao, XG Kong, Unique sequence 
characteristics of genes in the leftmost region of unique long region in 
duck enteritis virus. Intervirology. 52, 291-300 (2009). doi:1 0.1 159/000235742 
YF Li, B Huang, XL Ma, J Wu, F Li, W Ai, MX Song, HC Yang, Molecular 
characterization of the genome of duck enteritis virus. Virology. 391, 
151-152 (2009). doi:10.1016/j.virol.2009.06.018 

Y Zhao, JW Wang, F Liu, B Ma, Molecular analysis of US10, S3, and US2 in 
duck enteritis virus. Virus Genes. 38, 243-248 (2009). doi:1 0.1 007/s1 1262- 
008-0315-0 

Y Zhao, JW Wang, B Ma, FY Liu, Molecular analysis of duck enteritis virus 
US3, US4 and US5 gene. Virus Genes. 38, 289-294 (2009). doi:1 0.1 007/ 
S11262-008-0326-X 

Y Hu, H Zhou, Z Yu, H Chen, M Jin, Characterization of the genes encoding 
complete US10, SORF3, and US2 proteins from duck enteritis virus. Virus 
Genes. 38, 295-301 (2009). doi:1 0.1 007/s1 1262-009-0329-2 

Y Zhao, JW Wang, Characterization of duck enteritis virus US6, US7 and US8 
gene. Intervirology. 53, 141-145 (2010). doi:1 0.1 159/000274974 

A Dolan, FE Jamieson, C Cunningham, BC Barnett, DJ McGeoch, The 
genome sequence of herpes simplex virus type 2. J Virol. 72, 2010-2021 
(1998) 

DJ McGeoch, A Donald, FJ Rixon, Sequence determination and genetic 
content of the short unique region in the genome of herpes simplex virus 
type 1. J Mol Biol. 181, 1-13 (1985). doi:1 0.1 01 6/0022-2836(85)90320-1 
MF Summers, Zinc fingers motif for single stranded mucleic acids? 
Investigations by nuclear magnetic resonance. J Cell Biochem. 45, 41-48 
(1991). doi:1 0.1 002/jcb.2404501 1 0 

VR Holden, RR Yaamanchili, RN Harty, DJ O' Callaghan, Identification and 
characterization of an equine herpesvirus 1 late gene encoding a potential 
zinc finger. Virology. 188, 704-713 (1992). doi:10.1016/0042-6822(92)90525-T 
AJ Davison, DJ McGeoch, Evolutionary comparisons of the S segments in 
the genomes of herpes simplex virus type 1 and varicella-zoster virus. J Gen 
Virol. 67, 597-611 (1986). doi:1 0.1 099/0022-1 31 7-67-4-597 
BG Klupp, CJ Hengartner, TC Mettenleiter, LW Enquist, Complete, annotated 
sequence of the Pseudorabies virus genome. J Virol. 78, 424-440 (2003) 
EAR Telford, MS Watson, K McBride, AJ Davison, The DNA sequence of 
Equine herpesvirus-1. Virology. 189, 304-316 (1992). doi:1 0.1 01 6/0042-6822 
(92)90706-U 

H Delius, JB Clements, A partial denaturation map of herpes simplex virus 

type 1 DNA: evidence for inversions of the unique DNA regions. J Gen 

Virol. 33, 125-133 (1976). doi:1 0.1 099/0022-1 31 7-33-1 -125 

GS Hayward, RJ Jacob, SC Wadsworth, Anatomy of herpes simplex virus 

DNA: evidence for four populations of molecules that differ in the relative 

orientations of their long and short components. Proc Natl Acad Sci USA. 

72, 4243-4247 (1975). doi:1 0.1 073/pnas.72.1 1.4243 

K Hirai, K Ikuta, S Kato, Structural changes of the DNA of Marek's disease 

virus during serial passage in culture cells. Virology. 115, 385-389 (1981). 

doi:1 0.1 01 6/0042-6822(81 )901 1 9-7 

JW Balliet, PA Schaffer, Point mutation in herpes simplex virus type 1 oriL, 
but not in oriS, reduce pathogenesis during acute infection of mice and 
impair reactivation from latency. J Virol. 80, 440-450 (2006). doi:1 0.1 1 28/ 
JVI.80.1. 440-450.2006 



30. 



32. 



33. 



PC Chang, ML Hsieh, JH Shien, DA Graham, MS Lee, HK Shieh, Complete 
nucleotide sequence of avian paramyxovirus type 6 isolated from ducks. J 
Gen Virol. 82, 2157-2168 (2001) 

JD Parker, PS Rabinovitch, GC Burmer, Targeted gene walking polymerase 
chain reaction. Nucleic Acids Res. 19, 3055-3060 (1991). doi:1 0.1 093/nar/ 
19.11.3055 

M Reese, NL Harris, FH Eeckman, Large scale sequencing specific neural 
networks for promoter and splice site recognition, in Biocomputing 
Proceedings of the 1996 Pacific Symposium: 2-7 Jan 1996, ed. by Hunter L, 
Klein TE (Singapore: World Scientific Publishing Inc, 1996) 
SF Altschul, W Gish, W Miller, EW Mysers, DJ Lipman, Basic local alignment 
search tool. J Mol Bio. 215, 3389-3402 (1990) 



doi:1 0.1 1 86/1 743-422X-8-200 

Cite this article as: Liu et al:. Different linkages in the long and short 
regions of the genomes of duck enteritis virus Clone-03 and VAC 
Strains. Virology Journal 2011 8:200. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



